Collectivity-assisted ground state cooling of a nanomechanical resonator 
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We discuss cooling of a nanomechanical resonator to its mechanical ground state by coupling it 
to a collective system of two interacting flux qubits. We find that the collectivity crucially improves 
cooling by two mechanisms. First, cooling transitions proceed via sub-radiant Dicke states, and the 
reduced line width of these sub-radiant states suppresses both the scattering and the environmental 
contribution to the final phonon number. Second, detrimental carrier excitations without change in 
the motion of the resonator are suppressed by collective energy shifts. 

PACS numbers: 85.85.+j,85.25.-j,42.50.Wk,42.50.Nn 



Micro- and nanomechanical resonators (NAMRs) are 
currently in the focus of interest because of the fascinat- 
ing applications they promise [![. Access to these ap- 
plications requires cooling of the mechanical motion to 
the quantum mechanical ground state, which as yet is an 
unsolved task. Different routes towards the ground state 
have been suggested. For example, backaction cooling 
or sideband cooling in optomechanical systems has been 
proposed and demonstrated via cavity-assisted radiation 
pressure I n the electromechanical domain, cool- 

ing methods have been suggested which are based on 
auxiliary devices such as superconducting transmission 
lines 0], superconduc ting single-electron transistors Q, 
or quantum dots 0, Il0l | Just as in laser cooling of 
atoms or ions, backaction cooling or sideband cooling 
can only reach the ground state if the mechanical res- 
onance frequency exceeds the bandwidth of the cavity 
or auxiliary quantum system. This so-called strong con- 
finement regime (SCR) restricts the accessible parameter 
range, and is often difficult to implement. Active feed- 
back cooling [1, [n], [13] is a candidate for overcoming this 
limit. However, it typically requires difficult control and 
precise measurement of displacement of the resonator. 
Another route is to make use of quantum interference. 
In [l3| , a system was proposed in which the mechani- 
cal displacement of the oscillator changes the damping, 
rather than the frequency of a coupled cavity. This leads 
to destructive interference in the noise contributions, and 
thus to the possibility of ground state cooling already 
outside the SCR. Recently, we analyzed a cooling scheme 
for NAMR based on electromagnetically induced trans- 
parency (EFT) [bl ] . EIT cooling has originally been pro- 
posed [15( and demonstrated [161 ] for trapped ions, and 
works by eliminating unwanted heating transitions via 
destructive quantum interference. 

In this Letter, we propose an efficient ground state 
cooling scheme for a NAMR without counterpart in the 
cooling of atoms or ions. The NAMR is coupled to two 
interacting flux qubits as shown in Fig. [TJ The mutual 
interaction of the qubits gives rise to frequency-shifted 
collective qubit states with modified decay rates. We 
find that in a suitable cooling field configuration, the col- 



lective frequency shift effectively eliminates detrimental 
qubit carrier excitations without change in the NAMR 
motion. At the same time, cooling via the sub-radiant 
collective qubit state suppresses both contributions from 
the thermal environment and the scattering of the cooling 
field to the cooling limit. Due to these collective effects, 
our scheme offers efficient ground state cooling already 
for rather small driving fields, and for NAMR operating 
outside the strong-confinement regime. An implementa- 
tion is further assisted by the simple level structure of 
the qubit part compared, e.g., to EIT cooling, as the flux 
qubits arc modelled as two-level systems operating close 
to the optimum point. Unlike in backaction cooling and 
feedback cooling, the final NAMR state has no coherent 
shift in its phonon number. 

Before we proceed with the analysis, we start by ex- 
plaining the physical mechanisms behind collectivity- 
assisted cooling. First, we consider the two interact- 
ing qubits without coupling to the motional state |n) 
of the oscillator 17[. Qubit j (j G {L,R}) can be 



visualized as a two-level quantum system with ground 
state |<7j), excited state |ey), and decay rate 7. The 
two qubits L and R are coupled via their mutual in- 
ductances, such that the system dynamics can conve- 
niently be described in the collective-state basis \g) = 
\9l9r), \s) = (\e L gR) + |.9Lefl))/V2, |e) = \e L e R ), and 
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FIG. 1. (Color online) A nanomechanical resonator coupled 
to two flux qubits. The flux qubits interact with each other 
through their mutual inductance M and are damped via a 
common bath modeled as a LC circuit. 
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FIG. 2. (Color online) Energy level diagram for the NAMR 
coupled to two interacting flux qubits. Cooling (heating) 
transitions \n) — > \n ± 1) proceed via the narrow collective 
sub-radiant state \a) and are denoted by blue (red) arrows. 
Carrier excitation \n) — » |n) (purple arrow) is possible only 
via the strongly detuned symmetric collective state. 



I a ) = ( |eM7B.) - \gLe.R))/y/2. In a strongly interacting 
system lla,[l9|, the anti-symmetric state |a) is character- 
ized by a significantly reduced decay rate r a <C 7, while 
the symmetric state has an enhanced rate V s sa 27. At 
the same time, the two states arc shifted in energy due 
to the qubit-qubit interactions. Analogously, a driving 
field which is applied symmetrically to the two qubits will 
couple to transitions \g,n) <-> |s,n) and \s,n) «-> |e, n), 
but not to those via \a,n). Next, we consider in addi- 
tion the coupling to the motional eigenstates \n) of the 
NAMR. A relevant part of the total energy level spec- 
trum is shown in Fig. [2l It turns out that the NAMR 
motion couples anti-symmetrically to the two qubits, be- 
cause a movement of the oscillator always increases the 
loop area of one of the qubits, while decreasing the area of 
the other. Therefore, the symmetrically applied driving 
field still couples the ground state \g, n) to the symmetric 
state |s,n) without change in the motion, but it couples 
I g, n) to the anti-symmetric states |a, retl) if the motional 
state increases or decreases. The three leading excitation 
channels from the ground state are depicted as arrows in 
Fig. [2] for a field in resonance with the red sideband to 
state \a,n— 1). Then, as in sideband cooling @, this 
leads to a cooling of the NAMR, but with substantially 
inproved performance due to the scattering via collective 
states. The cooling limit crucially depends on the ratio 
of the NAMR eigenfrcquency v to the width of the scat- 
tering state. Sideband cooling to the motional ground 
state is only possible in the resolved regime in which 
v exceeds the scattering width. This explains the first 
advantage of collectivity-assisted cooling: The cooling 
limit is determined by Y a /v rather than 7/1/ as in most 
other cooling schemes. In particular, in the common case 
r a <C v < 7, a single qubit would be in the non-resolved 
regime, whereas the collective system effectively becomes 
resolved. The second advantage is a suppression of car- 
rier excitations \n) — ► which also leads to a heating. 
This suppression is achieved mainly by the energy shift 
of the collective states |s) and |a), which moves the sym- 
metric carrier scattering channel out of resonance with 



the cooling field. Thirdly, the narrow width r a also sup- 
presses the non- resonant heating processes |n) — > |n + l) 
as compared to the single-qubit case. These advantages 
together lead to a greatly improved collectivity-assisted 
cooling performance. 

We now proceed with a quantitative analysis. The 
NAMR has effective mass M c g, length I, frequency v 1 
and quality factor Q. It can be treated as a harmonic 
quantum resonator with Hamiltonian H r = hvb%. The 
quantized displacement operator is x = Xo(b + b^) 1 where 
b is the annihilation operator. Xo = ^Jh/2AI c gv denotes 
the zero point fluctuation of the NAMR. Each flux qubit 
consists of superconducting loops with three Josephson 
junctions, two identical and one junction smaller by a 
factor a. The qubits are exposed both to a constant 
magnetic field B perpendicular to the plane and to a 
common driving microwave field ^(i), which induces 
a time-dependent magnetic flux (TDMF). Following a 
standard treatment of the coupled flux qubits [2(| , both 
qubits are modeled as two-level systems. The free sys- 
tem Hamiltonian of flux qubits and NAMR is given by 
H a = hvtfb+hA(R ee -R gg ) + hA(R ss -R aa ). The qubit 
operators are defined as Rjk = \j)(k\ (j,k g {e, s,a, g}). 
The detuning A = luq — col is the difference of the qubit 
transition frequency ujq and the TDMF frequency ljl- A 
includes the always-on coupling and the collective shifts 
induced by the coupling of the two qubits 2(| 21 1. The 
qubits operate near their optimal point / = 0.5, and are 
assumed degenerate, which is reasonable if their transi- 
tion energy difference is much smaller than A. We assume 
that the qubit transition frequencies are large enough to 
neglect thermal excitations. 

The qubit-NAMR interaction is modelled by Hi = 

\/2M7(i? es +R S g +R.se +Rgs 

Rea — Rae)- The Rabi frequencies are Mlj 
with coupling energy aEj, KIl = KIr = 
(e^l sin(2</3p ' + 2Ttf)\gj) and Cj = (ej\cos(2tp 

„Cj) 
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2irf)\gj). The phase ip { p j) = (^ 
tribution from the momentum can be neglected. The 
Lamb-Dicke parameters evaluate to r/j = ^BIXqI-k j '$o> 
where <^ = Cj/Sj. We assume that 77 = t]l = t/r. 
Our system couples to an environment of temperature 
T, such that the thermal NAMR occupation is Ni = 
[exp^is/hsT) — 1] . In Born-Markov, rotating wave, 
and Lamb-Dicke approximations, the master equation is 



Sj = 



^)/2. 



p = - -[H + H u p] + % SEP + se^p 

+ [Ni + 1]JS? (y/Q, b) + Ni£> (v/Q, 6 f ) , (la) 

= J2 & ( r 0A \ej) {ej I - \gj){gj\) p, (lb) 
i 

■%?SEP =S£ (r s , R se + R gs ) + J£{T a , Rae — Rga) 

+J? (T s r) 2 ,(R ae -R ga )(b + tf)) 
+J? (T a r] 2 ,(R se + R gs )(b + tf)) 
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FIG. 3. (Color online) Final average phonon number n BB and 
its contributions as a function of the cooling field Rabi fre- 
quency fl. Parameters are Ni = 400, Q = 10 , F a = 0.l7,and 
i) = 3x 10~ 3 . Dephasing is (a) r = and (b) T = O.57. A 
fit to the environmental contribution by Eq. (J5| is shown as 
dashed line, with (a) C = 0.16 and (b) C = 2. 



+J8f (r s77 , (R ga - R ae )(b + tf),R se + R gs ) 
+J?{T aV , {R se + R gs )(b + fet), R ga - R ae ) , (lc) 

JSf(«, A)p = I {2ApA* - [A* A, p] + } , (Id) 
Jgf(K,Ai,A a )/> = f E { 2A ^ A 1 - . ( lc ) 

with the damping coefficients T s = 7 + 712 and T a = 
7 — 712 and k = i^/Q. Here, 7 and 712 are the Einstein A 
coefficient and the dipole-dipole cross damping rate, re- 
spectively, is the pure dephasing of a single qubit. In 
Eqs. ((I]), Jz?seP describes the spontaneous emission of the 
flux qubits. The additional pure dephasing in described 
by j£f0/9. The last two terms in Eq. ([Taj) result from the 
coupling of NAMR and thermal environment. 

As we start cooling from thermal equilibrium giving 
an initial phonon number Ni = N(v), the final phonon 
number n ss is determined by two contributions: one pro- 
portional to the initial phonon number and the other 
from the scattering of the cooling field, 



C(77,0,r a )r a zy 



N, 



-G(Ti,n,r a ) 



4;/ 



(2) 



Here we introduce coefficients C(r], 17, T a ) and G(r], fi, r o ) 
to correct the deviation of our approximate analytical 
expression from the numerical results. The analytical 
parts of Eq. ((2]) allow for a comparison with other cool- 
ing methods. We find that the environmental contribu- 
tion to the cooling limit is suppressed by a factor of CY a v 
in comparison with sideband cooling [9( and backaction 
cooling For this comparison, we assume effective 
Lamb-Dicke parameters defined as t\ c r = rjfl/j in our 
case and rj e ff = r)y/n max in the backaction cooling scheme, 
where rj = (dui c / dx)(Xo/i>) with the cavity resonance fre- 
quency lo c and n max as the photon number in resonance 
[6J. Interestingly, as the environmental contribution is 
proportional to is, the nonresolved regime with small v 
is favourable. We also find that in our cooling scheme 
there is no coherent shift in the final phonon number, 
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FIG. 4. (Color online) The final resonator phonon number n ss 
[curves (i) and (iv)] and its contributions from environment 
(ii) and scattering (iii). The parameters are = and 
SI = 47 in curves (i)-(iii), and = O.57 and = IO7 in 
(iv). Further, r/ = 3 x 10" 3 , and Q = 10 6 . Dashed lines 
show analytical fits to the numerical data using Eq. ([2]). (a) 
shows the results versus the sub-radiant decay width F a with 
v — O.57, and (b) against the resonator frequency v with 

r a = 0.057. 



other than in backaction cooling [6| and feedback cool- 
ing [T2I ]. This is crucial for many applications such as in 
quantum information science [221 ] - 

We now turn to a full numerical study of Eqs. UJ. 
Typical flux qubit parameters we envisage are a = 
0.8, Ej/Ec = 100, Ej = 1000 GHz (I c - 300 nA), where 
Ej and Ec are the Josephson and charging energies of 
junction, respectively. We assume a bias flux /& = 0.5005. 
We calculate transition frequency ujq/Ej ~ 4.5 x 10~ 3 , 



and transition matrix element 5, 



-0.16 and ? - 0.05. 



We assume a decay rate of about 2n x 2 MHz, and ex- 
perimental results show that this rate can remain similar 
over a wide range of bias [23] . Pure dephasing is assumed 
as half of the decay rate for a t ypi cal noise spectral den- 



sity (/i$ ) /Hz at f b = 0.5005 [23j. Note that the qubit 



parameters can be engineered over a wide range [24| . Re- 
garding the resonator, a double clamped NAMR with size 
of 300 /imxlOO nmxlOO nm made from silicon nitride 
can be manufactured with frequency v ps 27t x 1 MHz 
and quality Q w 10 6 j25|. This NAMR leads to a Lamb- 
Dicke parameter r\ ~ 0.003 for B < 2T. Thus our system 
can be realized based on NbN superconducting qubits. 
Our sample chip could be placed in a dilution refriger- 
ator of 20 mK giving rise to a initial phonon number 
Ni ~ 400. Assuming that the qubits is driven by a 10011 
microwave line coupled by mutual inductance ~ 50 fH, 
a Rabi frequency of IO7 only requires ~ 0.5 nW input 
power. This small input power helps to reduce nonlin- 
ear and phase noise [7|. Throughout our investigation, 
we fix A = 5OO7 and apply the TDMF with detuning 
A = v + A corresponding to a red sideband excitation 
\g,n) — ► I a, n— 1). Note that a large mutual inductance 
M can lead to couplings A up to a few GHz 



17] 



An example for the dependence of the final phonon 
number n ss on the TDMF strength fl is shown in Fig. [3] 
Without pure dephasing, one can ideally cool a NAMR 
from Ni = 400 to a final phonon number n ss smaller than 
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0.05 in the range of 3 < Sl/7 < 8, see Fig.[3]Ja). This cor- 
responds to a ratio Ni/Nf of initial to final phonon num- 
ber exceeding 10 3 . By fitting the environmental contribu- 
tion, we find that C = 0.16. Our analytic formula Eq. @ 
becomes invalid for O exceeding 47. Taking into account 
an additional pure dephasing = O.57, we find that 
n ss is still lower than 0.5 in the range of 4 < SI/7 < 20, 
see Fig. EJb). At SI = IO7, the ground state occupation 
is 76%. Further investigation show that a ground state 
cooling is possible as long as < 7. 

According to Eq. @, both the environmental and 
the scattering contributions to the cooling limit are sup- 
pressed by the reduced decay rate r a of the sub-radiant 
state. This result is confirmed by our numerical results 
shown in Fig.JUa). Interestingly, an optimal cooling limit 
is obtained for r a w O.O57. In this point, the total final 
phonon number is only 0.03, which corresponds to a 97% 
occupation of ground state. Further decrease of T a leads 
to an unexpected rapid increase of the cooling limit. The 
origin of this increase is that due to level shifts, the qubit 
transition frequency is slightly off-resonant with the driv- 
ing field frequency. If the scattering state becomes too 
narrow, then the driving field cannot excite the cooling 
transitions any more. We verified this interpretation by a 
suitable slight change in the driving field frequency, which 
allows to achieve effective cooling also for smaller T a . 
With additional pure dephasing = O.57, the optimum 
cooling limit is obtained for SI = IO7. For r a = 0.057, 
the steady state is n ss w 0.29, which indicates 78% occu- 
pation of ground state. 

Next, we study the dependence of the cooling limit on 
the resonance frequency v of the NAMR, sec Fig. 0|b). 
The environmental contribution to n ss is proportional 
to v. Our scheme thus provides particularly efficient 
cooling of NAMRs with low frequencies. But since 
the scattering contribution is inversely proportional to 
the square of v, this part becomes eventually domi- 
nant with decreasing v. We find an optimal point of 
v ~ 0.257 where a minimum final phonon number of 
0.02 is achieved. This number agrees well with the op- 
timal value v opt = \J rj 2 Q, 2 QT a G/8CNi evaluated from 
Eq. (0) with C = 0.4 and G ~ 3. When including a 
pure dephasing of — O.57, an optimal cooling limit of 
n ss = 0.15 is obtained for SI = IO7 and v = 7. 

In summary, we have presented an efficient ground 
state cooling scheme for a NAMR which operates by cou- 
pling the NAMR to two interacting flux qubits. We found 
that collectivity enhances the cooling performance in two 
crucial ways. First, the collective frequency shift of the 



qubit states effectively eliminates unwanted carrier exci- 
tations. Second, cooling via the sub-radiant Dicke state 
suppresses both contributions from the thermal environ- 
ment and the scattering of the cooling field to the cooling 
limit. As a result of this, ground state cooling can be 
achieved already in the non-resolved regime, and using 
rather small driving fields. 
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